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I.  INTRODUCTION 

Pattern  synthesis  methods  exist  for  many  different  antenna  arrays 
under  various  constraints  [1].  In  this  report  a radiation  field  pattern 
magnitude  synthesis  procedure  for  an  N-port  reactlvely  loaded  circular 
dipole  antenna  array  Is  presented.  The  procedure  applies  an  optimiza- 
tion algorithm  directly  to  the  synthesis  error  function  to  find  the 
reactive  loads  required  for  minimum  error.  It  accepts  any  arbitrary 
magnitude  pattern  as  an  Input. 

In  many  physical  cases  only  the  far  field  power  pattern  Is  speci- 
fied. This,  In  effect,  specifies  the  far  field  magnitude  and  leaves  the 
far  field  phase  arbitrary.  Normally,  the  far  field  phase  does  not  affect 
system  performance.  The  synthesis  procedure  used  In  this  report  requires 
a specification  of  the  magnitude  of  the  far  field  electric  field  with  no 
restrictions  on  the  phase.  Since  the  far  electric  field  for  the  circular 
array  considered  in  this  report  Is  linearly  polarized,  the  square  of  the 
far  field  electric  field  Is  proportional  to  the  far  field  power  pattern. 

The  radiation  pattern  of  an  antenna  array  can  be  controlled  by  placing 
reactive  loads  at  the  antenna  ports.  The  reactlvely  loaded  circular  di- 
pole array  to  be  considered  has  the  center  dipole  fed  directly  and  the 
outer  dipoles  parasltically  excited  with  reactive  loads  across  their  In- 
put terminals  [2].  The  advantage  of  reactlvely  loaded  antenna  arrays  Is 
the  elimination  of  a transmission  line  feed  system  to  all  of  the  elements. 
The  circular  array  geometry  was  chosen  because  of  Its  N-fold  rotational 
symmetry  with  respect  to  4>  In  a constant  0 plane  (see  Figure  1).  There- 
fore, an  arbitrary  synthesized  pattern  In  a constant  6 plane  is  easier 
to  realize.  A reactlvely  loaded  linear  array  has  much  less  symmetry 
which  makes  it  harder  to  synthesize  an  arbitrary  pattern. 
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II.  GENERAL  FORMULATION 

Consider  a source  1 and  the  field  f it  produces  on  the  radiation 
sphere.  Also,  let  T be  a knovm  operator  representing  the  antenna  system. 

The  analysis  problem  is  concerned  with  determination  of  the  radia- 
tion characteristics  f for  a given  current  distribution  i of  an  antenna 
T,  that  is. 


Ti  = f (1) 

where  this  is  an  exact  relationship. 

The  synthesis  problem  is  concerned  with  determining  the  current  dis- 
tribution i of  an  antenna  T given  a specified  field  pattern  f^,  that  is, 

Ti  ==  f . (2) 

o 

This  is  usually  an  Inexact  relationship.  Deschamps  and  Cabayan  [3]  dis- 
cuss the  general  Inverse  problem  and  the  difficulties  encountered. 

In  order  to  use  matrix  operations  in  the  synthesis  problem,  (2)  is 
discretized  as  follows.  Let 


i = I i e (3) 

n n 

n=l 

where  i are  constants  and  e are  basis  elements.  Since  the  source  is 
n n 

continuous,  the  e are  functions  and  (3)  is  an  approximation  to  the  true 
source.  Define  i to  be  the  vector  having  the  components  i^,  that  is, 

^ = f^n^Nxr 

Next,  substitute  (3)  into  (2)  and  evaluate  (2)  at  P points  (0p»<i>p).  P = 
1,  2,  ...,  P on  the  radiation  sphere.  In  matrix  notation,  we  have 

[T]T  ^ f (5) 


where  ? is  the  vector 
o 


fo=  ffopl 
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and  [T]  is  the  matrix 

tT]  = KTe^)p]p,N  • (7) 

Here  f denotes  the  value  of  f at  the  point  (0  ,(p  ) and  (Te  ) denotes 
op  o p p n p 

the  pattern  of  e evaluated  at  the  point  (6  ,(p  ) . 

n p p 

Next  consider  an  array  of  N z-directed  center-driven  parallel  one- 
half  wavelength  long  dipoles.  The  boundary  condition  to  be  satisfied  on 
the  surface  of  the  mth  dipole  is 

-E^(z)  = V 6(z)  (8) 

m m 


where  E (z)  is  the  z component  of  the  tangential  electric  field  induced 
m 

by  the  electric  currents  on  the  dipoles,  V is  the  voltage  applied  at 

m 

the  center  of  the  mth  wire,  and  6(z)  is  the  Dirac  delta  function.  The 
electric  current  on  the  nth  wire  is  approximated  by  a filament  current 

I cos  kz  (9) 

n 

where  I Is  a constant  to  be  determined  and  k is  the  propagation  constant, 
n 

Substituting  (9)  into  (8),  we  obtain 

N 

- y IE  (z)  = V 6(z)  (10) 

n mn  m 

n=l 

where  I E (z)  is  the  z component  of  the  electric  field  induced  by  (9)  on 
n mn 

the  surface  of  the  mth  dipole.  Following  Galerkin's  method,  we  multiply 
(10)  by  cos  kz  and  Integrate  from  -X/4  to  X/4  (X  is  the  wavelength)  on  z 
to  obtain  the  matrix  equation 

[Z]f  = V.  (11) 


Here  [Z]  is  the  impedance  matrix  of  the  antenna  array  with  elements 

X/4 
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For  the  reactlvely  loaded  circular  dipole  array,  (11)  becomes 


where 


[Za  + ZJI  = V 


[Za]  = [Z] 


(13) 


(U) 


Is  the  port  Impedance  matrix  for  the  circular  dipole  array  with  the 
elements  defined  by  (12), 

0 0 . . . 0 


0 


0 


[o  0 . . . jX^ 


(15) 


Is  the  reactive  load  diagonal  matrix  where  jX^^  = 0 corresponding  to  the 
central  dipole. 


V = 


(16) 


is  the  excitation  voltage  vector  with  one  non-zero  element  V corresponding 
to  the  central  dlpole» 


I 

L n 


(17) 


is  a column  vector  of  the  port  currents. 


1 


The  far  electric  field  for  the  N dipoles  in  the  (p  direction  and 
0 = tt/2  plane  is 

-jkr  N jk(x  cos(fi+y  sin(J)) 


where  n is  the  intrinsic  Impedance  of  free  space, 


-in  e 


2it  r 

N jk(x  cos<t>+y  sin(})) 

FOP)  = I I e " 


Since  we  are  Interested  only  in  the  pattern  shape,  the  far  electric  field 

E (6)  can  be  represented  by  F(<p). 
z 

The  pattern  magnitude  synthesis  procedure  uses  P points  from  the 
desired  pattern  |f  | and  finds  a source  I such  that  the  pattern  error 

P P N 

C - I w l|F«  )|  - )||'  ■ L W II  I I - |£„»  )||^  (21) 

p=l  p=l  n=l 


is  minimized.  Here  w^  is  a weighting  coefficient  (w^  > 0) , 


T = e 
pn 


ik(x  cos(p  +y  sinif  ) 
n P u P 


I = {[Z.  + Z,  ]~  V}  . 
n A L n 


The  { ) in  (23)  indicates  the  nth  element  of  the  column  matrix 
n ^ 

[Z  + Z ] V.  Substituting  (23)  into  (21),  we  obtain 
A L 


G = I w II  I {[Z^  + ZJ~  V)  T I - |f„(^  )il  . 


p=l  ^ n=l 


There  are  N unknowns  in  (24),  namely,  N-1  unknown  reactive  loads 
{X^}  and  one  unknown  excitation  coefficient  V^.  If  [Z  ] could  be  separated 


i 


7 

out  from  the  Inverse  operator  in  (24),  an  exact  expression  could  be  ob- 
tained for  determining  the  required  loads  ix")  for  minimum  error.  Since 
this  Is  not  mathematically  feasible,  an  approximation  technique  is  used 
for  determining  the  {X^}  and  then  incrementing  these  values  in  a direction 
which  minimizes  (24).  One  pos.sibie  approach  is  the  following. 

Let 


\ = iX,  + (25) 

where  is  an  incremental  change  in  . Then  [Z^  + Z^  ] ^ can  be  expanded 
using  a Maclaurin  series  as 


[Xa  + jX|  + Zj]  ^ [I  + A + + ...][Za  + jX^]"^  (26) 


where 

A = -[Za  + jXj  )"^[Zj  J (27) 

I = Identity  matrix. 


Assuming  (z^  1 are  small  and  the  , A , ...  terms  in  (26)  are  negligible 
in  relation  to  [I  + A],  we  can  write  (26)  as 


[X-A  + iXj  + Zj  ] 1 = [I  - [Za  + jX^  ]][Za  = (28) 


Substituting  (28)  into  (24),  we  obtain 
P N 

1 I 


- I “II  I fit  - IZ,,  + JX.r'ujUZA  + 


pA  P n-1 


pn' 


Po<Vl'  • 

(29) 


If  the  , A , ...  terms  in  (26)  are  not  negligible  compared  to  [I  + A], 
new  fX^}  must  be  chosen  and  small  changes  made  in  z^  to  test  for  a minimum. 
(Note  that  since  (29)  contains  absolute  value  operations,  the  minimum  does 
not  necessarily  occur  at  a stationary  point.) 

The  advantage  of  this  technique  lies  in  the  fact  that  [z^]  in  (29) 
can  be  varied  without  performing  an  inverse  operation.  A disadvantage  is 
the  fact  that  the  {X^}  have  to  be  chosen  so  that  A^,  A^,  ...  are  suffi- 
ciently small  for  (28)  to  hold.  Also  lacking  is  a good  strategy  for  in- 
crementing the  fX?)  in  the  direction  of  a minimum. 
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Another  method  of  minimizing  (24)  is  the  direct  application  of  an 
optimization  algorithm.  This  approach  Is  used  in  this  report  since  the 
method  Is  not  dependent  on  the  proximity  of  the  to  a minimum  point, 

and  the  optimum  seeking  strategy  for  Incrementing  the  variables  Is  an 
Integral  part  of  the  algorithm. 

The  optimization  subroutine  selected  is  a Rosenbrock  algorithm  [4]. 
The  Rosenbrock  algorithm  uses  a set  of  N mutually  orthogonal  directions 
In  each  cycle  of  searches  (stage).  After  each  stage,  a new  set  of  ortho- 
gonal directions  are  determined  by  rotating  the  former  direction  vectors 
until  they  are  oriented  toward  the  direction  of  fastest  de  rease  of  the 
objective  function.  A description  of  the  Rosenbrock  algo  Ithm  along  with 
a computer  listing  of  the  subroutine  Is  included  in  the  Appendix. 

In  this  report  N variables  are  used  In  the  Rosenbrock  subroutine  — 
N-1  reactive  loads  and  one  excitation  voltage  coefficient.  It  should  be 
mentioned  that  since  (24)  does  not  depend  on  the  phase  of  V^,  can  be 
taken  as  real  and  positive.  When  this  is  done,  (24)  becomes  a simple 
quadratic  function  in  V . From  this  quadratic,  the  optimum  value  of  V 
can  be  easily  obtained  for  fixed  Z^.  If  this  optimum  value  of  V were 
substituted  back  into  (24),  only  the  N-1  reactive  loads  would  enter  Into 
the  Rosenbrock  search. 


III.  OPTIMALLY  EXCITED  ARRAY 

In  this  section,  a solution  is  given  for  the  same  array  analyzed  In 
the  preceding  section,  but  with  all  of  the  elements  driven  and  no  reactive 
loading.  In  the  results  section,  a compar i •■'^on  of  the  synthesized  patterns 
will  be  made  between  the  reactively  loaded  array  and  the  optimally  excited 
array  to  indicate  the  degree  of  degradation  due  to  the  restricted  excita- 
tion for  the  reactively  loaded  array. 

The  far  electric  field  magnitude  at  a point  p in  the  6 = tt/2  plane 
has  the  same  form  as  In  the  reactively  loaded  case,  that  is. 


F($p) 


N 


I 

n=  I 


IT  I . 

n pn  ' 


(30) 


Here,  i can  take  on  any  complex  value,  in  contrast  to  the  restricted 
n 


values  for  the  reactively  loaded  case.  A solution  for  the  pattern  magni- 
tude synthesis  problem  using  (30)  can  be  found  in  [5],  and  the  following 
is  a summary  of  the  formulation. 

The  error  function  to  be  minimized  for  the  general  synthesis  problem 
where  both  magnitude  and  phase  of  the  desired  pattern  ? are  specified  is 


P N 


1 w|2^  IT  -f((}))|. 

^ n ' ^ n ■nn  n n * 


Pi  n pn  op 
p=l  ‘ n=l  ‘ 


The  solution  to  (31)  (see  section  V of  [5]  with  a = 0)  is 


-t-  -1 

1 = [T  UT]  [T  W]f 


where  * signifies  complex  conjugate,  'signifies  transpose  and  [W]  Is  a 


weighting  matrix. 

When  the  pattern  magnitude  only  is  specified,  the  error  function 


becomes 


1 w iT  |-  f(())) 

p ' ' n pn ' 'op 

p=l  ^ n=l 


To  circumvent  the  inner  magnitude  operation  in  (33) , we  consider  the 
following  more  general  equation  which  has  a minimum  less  than  or  equal 
to  that  of  (33) , 


g(6)  = y w i y i T - If  (4)  )|e  p|‘ 


p=l  P n=l  P"  ° P 


In  other  words,  we  are  specifying  a phase  for  the  pattern  magnitude  in 


order  to  obtain  a more  easily  solved  expression.  For  i^  fixed,  a minimum 


is  obtained  when  both  terms  within  the  outer  magnitude  signs  of  (34)  are 


in  phase,  that  is. 


I i T 

- n^l  ^ P^ 


N 

II  i T I 

n=l  " P" 


I 


■? 

j- 


t.' 
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If  condition  (35)  Is  used  In  minimizing  (34),  the  minimum  of  (34) 
will  be  equal  to  the  minimum  of  (33).  (This  can  be  shown  by  substituting 
(35)  into  (34),  and  after  factorization  It  Is  evident  that  the  magni- 
tude of  the  resulting  expression  will  be  equal  to  that  of  (33).) 

To  find  the  minimum  of  (34),  the  following  Iterative  procedure  Is 

used : 

1.  Assume  starting  values  for  3.,  3« , ••.,  3 • 

^ P 

2.  Keep  the  B fixed  and  calculate  the  1 which  minimizes 

P 1 n 

E using  (32) . 

3.  Keep  the  1 fixed  and  calculate  the  3 which  minimizes 

n P 

E using  (35). 

4 . Go  to  step  (2) . 

This  procedure  converges  to  a stationary  point  because  steps  (2)  and  (3) 
cannot  Increase  e. 


IV.  NORMALIZED  SYNTHESIS  ERROR  AND  Q-FACTOR 

Two  figures  of  merit  will  be  used  in  the  next  section  to  evaluate 
the  synthesis  results.  The  first  figure  of  merit  is  the  normalized  synthesis 
error  which  Is  defined  as 


e 

syn 


y w \ \¥(<p  ) I - |f  ((^  ) 

D P O p 


P=i  P 


y w if  ((j) ) 

p " p 

p-i 


(36) 


^It  should  be  noted  that  if  [f*WT]  is  ill-conditioned,  the  procedure 
will  fail  since  small  rounding  errors  that  occur  In  its  inversion 
can  cause  large  errors  in  its  inverse  [T*WT]  To  use  this  pro- 

cedure, the  ratio  of  the  magnitude  of  the  largest  eigenvalue  to  the 
smallest  of  [T*WT]  should  be  on  the  order  of  10  or  less. 


I 
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The  smaller  the  magnitude  of  the  closer  the  synthesized  pattern 

Is  to  the  desired  pattern. 

The  second  figure  of  merit  Is  the  Q or  quality  factor.  The  Q- 
factor  Is  a measure  of  the  magnitudes  of  the  element  excitation  required 
to  produce  a given  pattern.  The  larger  the  value  of  Q,  the  more  Imprac- 
tical the  pattern  realization  becomes  due  to  the  large  element  excita- 
tions required.  Lo,  Lee,  and  Lee  defined  the  Q-factor  as  [6]: 


I |J  I 

n=l 

|p(0,4))l^dfi 


where  J Is  the  nth  element  current  excitation  value  and  P(0,(J))  Is  a 
n 

pattern  function, 

A numerical  approximation  to  (37)  Is 


Q = P 


I ) 

p=i 


The  multiplier  P is  Introduced  to  make  Q relatively  insensitive  to  the 
number  of  field  points  chosen.  The  main  difference  between  the  Q ex- 
pressions of  (37)  and  (38)  is  that  the  denominator  of  (37)  Is  an  Inte- 
gral over  the  whole  radiation  sphere  while  the  denominator  of  (38)  Is 
a summation  over  P points  In  the  XY  plane. 


V.  klfreslntative  computations 

A computer  program  has  been  written  using  the  equations  derived 
in  the  preceding  sections.  A description  and  listing  of  the  optimiza- 
tion subroutine  is  included  in  the  Appendix  of  this  report.  In  this 
section  results  are  given  to  illustrate  the  synthesis  technique  applied 
to  a reactively  loaded  array  for  various  pattern  shapes.  For  all  cases 
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except  the  reduction  of  an  unwanted  side  lobe  (Figure  13)  all  w 

P 

used  are  equal  to  one. 

Figure  2 Illustrates  the  antenna  configurations  for  both  the 
seven  element  (array  A)  and  the  thirteen  element  (array  B)  circular 
arrays  to  be  used  in  the  synthesis  procedure. 

Figures  3 to  6 Illustrate  Gaussian  shaped  synthesis  patterns  in 
the  (J)  = 0 and  (})  = 30°  directions  for  reactlvely  loaded  arrays  A and  B 
In  which  the  specified  pattern  beam  widths  are  varied.  For  all  cases 
array  B (with  a greater  number  of  elements  than  array  A)  realized 
lower  synthesis  errors  than  array  A.  The  worst  case  (largest  synthe- 
sis error)  of  the  three  pattern  shapes  to  be  synthesized  for  both 
arrays  A and  B is  pattern  (a).  Pattern  (a)  Is  the  narrowest  of  the 
three  patterns  and  apparently  a greater  number  of  elements  than  thir- 
teen Is  needed  to  realize  the  directivity  required  by  this  pattern 
shape.  Pattern  (b)  is  synthesized  with  a slightly  higher  synthesis 
error  than  pattern  (c)  for  array  A but  the  reverse  is  true  for  array  B. 
Overall,  pattern  (b)  appears  to  be  the  best  compromise  In  terms  of 
pattern  synthesis  error  and  distortion  when  specifying  a pattern  shape 
to  be  synthesized  for  arrays  A and  B. 

Tables  1-4  list  the  reactive  loads  required  for  the  patterns  shown 
In  Figures  3-6.  Note  that  loads  which  affect  the  beam  width  In  one 
direction  (therefore  changing  their  values  for  each  beam  width  varia- 
tion) do  not  necessarily  affect  the  pattern  In  another  direction.  For 
instance,  the  reactive  load  for  element  four  In  Table  1 varies  from 
+26  ohms  to  +90  ohms  as  the  pattern  beam  width  changes  In  the  4>  = 0 
direction  while  the  reactive  load  for  element  four  In  Table  3 does 
not  change  at  all  as  the  pattern  beam  width  changes  in  the  ({>  = 30° 
direction . 

Figure  5 illustrates  an  important  property  (necessary  for  syn- 
thesizing patterns  in  arbitrary  directions)  of  the  circular  array  due 
to  its  symmetry,  that  is,  a pattern  in  the  <f>  = 0 direction  can  be 
shifted  exactly  (for  arrays  A and  B)  to  the  ([)  = n60° , n = 1,  2,  ...,  5 
directions  by  cyclically  changing  the  reactive  loads  required  for  the 
0=0  direction.  The  pattern  in  Figure  7 is  obtained  by  taking  the 
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(a)  Array  A - N = 7 
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/ \.'V  ;\ 

/ .A  '60°  \ 


(h)  Array  B - N = 13 

Fi^.  2.  Antenna  configuration  for  a seven  element  (array  A) 
and  thirteen  element  (array  B)  circular  array. 


Far  electric  field  pattern  variations  for  reactively  loaded  circular  array  A where 
A - pattern  shape  to  be  synthesized,  and  element  spacing  s = 0.25X. 


pattern  variations  for  reactively  loaded  circular  array  B where 
be  synthesized,  and  element  spacing  s = 0.25X. 
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reactive  loads  for  pattern  1 of  Figure  4 and  shifting  them  accordingly. 
For  Instance,  to  shift  the  pattern  from  the  (p  = 0 to  the  (J)  = 60°  direc- 
tion, the  reactive  load  value  for  element  two  becomes  the  reactive 
load  value  for  element  four;  the  reactive  load  value  for  element 
three  becomes  the  reactive  load  value  for  element  five,  etc. 

Figures  8 to  10  illustrate  synthesized  patterns  having  a single 
main  lobe  for  both  optimally  excited  (used  for  comparison)  and  reac- 
tlvely  loaded  arrays. 

In  Figure  8 a Gaussian  shaped  pattern  similar  In  beam  width  to 
that  of  pattern  (b)in  Figures  3 to  6 Is  synthesized.  Both  excitations 
achieve  a low  synthesis  error  with  about  the  same  Q. 

In  Figure  9,  a sector  pattern  is  synthesized.  Since  this  is  not 
a natural  pattern  shape,  each  excitation  type  attempts  to  distort  the 
normal  beam  shape  to  fit  the  contour.  The  synthesis  error  for  both 
excitation  types  is  still  low  due  to  the  relatively  large  beam  width 
required  (realizable  as  Illustrated  in  Figures  3 to  6). 

In  Figure  10,  a rectangular  shaped  pattern  is  synthesized  which 
requires  a narrow  beam  width.  The  synthesis  error  Is  higher  for  this 
case  than  In  the  previous  two  cases.  A linear  array  configuration 
would  probably  realize  a lower  synthesis  error  than  the  circular  array 
configuration  due  to  the  Increased  directivity  characteristics  in  the 
end  fire  direction. 

An  attempt  to  synthesize  multiple  main  lobes  is  shown  In  Figures  11 
and  12.  The  synthesized  patterns  obtained  for  array  B are  much  better 
than  for  array  A.  This  degree  of  improvement  is  not  evident  in  Figures 
8 to  10.  It  appears  that  more  elements  are  needed  to  synthesize 
multiple  shaped  beam  patterns. 

The  pattern  illustrated  in  Figure  13b  is  an  attempt  to  synthesize 
the  sector  shaped  pattern  without  the  side  lobe  indicated  by  the  arrow 
In  Figure  13a.  This  is  accomplished  by  weighting  the  synthesis  error 
greater  in  the  arrow  direction  than  in  all  other  directions.  As  shown 
In  Figure  13b,  the  unwanted  side  lobe  is  decreased  at  the  expense  of 
increased  side  lobes  in  other  directions  and  greater  overall  synthesis 


error . 


(b)  Array  B 

Far  electric  field  synthesized  patterns  for  a reactively 
loaded  circular  array  with  element  spacing  s = 0.35X 
(A-pattern  shape  to  be  synthesized,  OK-optimally  excited 
RI, -reactive  loading). 


(a)  Array  A 


Fig.  11.  Far  electric  field  synthesized  patterns  for  a reactively 
loaded  circular  array  with  element  spacing  s = 0.35X 
(/■-pattern  sliape  to  he  synthesized,  OE-optimaliy  excited 
RL-reactive  loading). 


(a)  Array  A 


(b)  Array  H 

f’ig.  12.  Far  electric  field  synthesi^.ed  patterns  fcir  a reactively 
loaded  circular  array  with  element  spacing  s = 0.35X 
(A-pattern  shape  to  be  synthesized,  OF.-opt imal ly  excited, 
RL-reactive  loading). 
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(a)  Synthesized  pattern  with  unwanted  side  lobe. 


(b)  Synthesized  pattern  with  unwanted  side  lobe  removed. 

Fig.  13.  Synthesized  pattern  with  side  lobe  restriction  for  array  B 

with  element  spacing  s = 0.35X  (A-pattern  shape  to  be  synthe 
sized,  RL-reactlve  loading,  \ -unwanted  side  lobe  direction) 


VI. 


CONCLUSIONS 


The  antenna  element  currents  for  the  reactively  loaded  array  are 
a subset  of  currents  for  the  optimally  excited  array.  Therefore,  while 
the  error  in  the  reactively  loaded  array  patterns  approaches  that  of  the 
optimally  excited  array,  it  is  always  greater.  The  patterns  are  closest 
when  the  specified  shape  is  typical  of  the  patterns  realizable  by  re- 
actively loaded  arrays  (pattern  (b)  in  Figures  3-6).  Verification  of 
this  statement  is  sliowm  in  Figure  8. 

The  synthesis  technique  used  in  this  report  is  adaptable  to  problems 
Involving  more  variables,  like  antenna  element  length  and  position.  Also, 
various  weighting  functions  can  be  used  to  improve  tiie  convergence  of  the 
optimization  algorithm. 

Future  work  should  investigate  more  fully  the  correlation  between 
antenna  array  configuration  and  synthesis  pattern  shapes  for  reactively 
loaded  arrays.  These  results  would  allow  specifying  an  antenna  array  con- 
figuration for  a corresponding  pattern  shape  before  the  synthesis  procedure 
has  begun . 
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APPENDIX 


A.  THE  SUBROUTINE  ROSIEM 
1)  Theory : 

The  Rosenbrock  search  technique  uses  N mutually  orthogonal 
directions  during  eacli  search  cycle  to  find  a relative  minimum.  This 
strategy  differs  from  a steepest  descent  technique  which  uses  successive 
orthogonal  directions,  but  these  successive  directions  do  not  necessarily 
form  a mutually  orthogonal  set. 

The  basic  elements  of  the  Rosenbrock  algorithm  are  as  follows: 

a)  Step  size  - The  step  size  in  a given  direction  is  chosen 

by  specifying  an  arbitrary  magnitude  e and  then,  if  a step 
e decreases  the  value  of  the  function  (for  a minimization 
problem),  e is  multiplied  by  a constant  a (a  > 1).  If  the 
value  of  the  function  increases,  e Is  multiplied  by  a con- 
stant -8  (0  < 8 1).  e,  a,  and  3 have  to  be  specified  by 

the  user.  e should  be  on  the  order  of  one  percent  of  the 
average  magnitude  of  the  variables.  To  determine  the  values 
of  a and  3,  a series  of  trials  using  various  values  of  n 

and  3 should  be  made  for  a given  class  of  functions  (preferably 
for  a known  function  for  which  a solution  is  available)  for 
wliich  tlie  function  to  be  minimized  belongs. 

b)  Direction  - The  Rosenbrock  algorit.im  uses  N orthogonal  direc- 
tions d , d , ....  d at  each  stage  rather  than  choosing  a 

12  N 

single  direction  in  whicli  to  progress.  A search  is  made  in 
each  orthogonal  direction  before  continuing  to  the  next  stage 
(at  least,  one  trial  has  been  successful  (a  value  less  than 
or  equal  to  tlie  old  value)  and  one  lias  failed  in  each  direction). 
There  are  three  cases  to  consider: 

i)  The  first  trial  is  a success, 

ii)  Tlie  first  trial  is  a failure  and  the  second  trial  is  a 
success. 


iii)  Tlie  first  trial  is  a failure  and  the  second  trial  is  a 
f ai lure . 
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For  a start,  let  [D^]  = [I]  (Identity  matrix)  where  d^,  d^, 
are  individual  directions  in  which  steps  are  to  be  taken.  The  search 
technique  progresses  as  follows  (for  one  variable  x,): 

1)  A step  is  taken  in  the  direction  d^  from  an  initial  point 

^ 0 ,0. 

2)  If  f(x^)  < f(x^),  the  step  Is  successful  (case  (1))  and  the 
step  size  Is  multiplied  by  a again  and  again  until  a failure  Is  recorded. 
After  this  failure,  the  result  of  the  last  successful  trial  and  the 
number  of  successful  trials  Is  stored. 

3)  If  f(x^)  > f(x*^)  on  the  first  trial  (a  failure  - case  (11)),  the 
sign  of  the  step  size  e is  changed  and  the  search  continues  as  In  (2). 

4)  If  f(x^)  > f (x^)  and  f(x^)  > f (x*^)  , the  second  and  further 
trials  are  determined  by  multiplying  e by  minus  6 until  a success  Is 
recorded . 

After  a search  has  been  made  for  each  variable  in  all  directions 
d^,  a new  [D^]  is  determined  using  the  Gram-Schraldt  orthogonallzatlon 
procedure.  The  orthogonallzatlon  procedure  completes  the  first  search 
cycle  and  the  same  procedure  Is  repeated  in  the  new  direction  [D^].  The 
procedure  is  continued  until  a convergence  criterion  Is  satisfied. 


2 ) Descript  i oti 

The  subroutine  ROSim  (N, N1,N2,.AA,BB, STEP ,CUL,  EF, ERR, GOM, 
PAT,V,XLOAD,ZA)  determines  the  loads  } and  excitation  voltage 
which  minimize  the  pattern  synthesis  error  e (equation  24). 

The  input  variables  are  defined  as  follows: 

N = number  of  antenna  elements 
N1  = number  of  search  cycle  (if  N1  = 1,  [D]  = [I]  ) 

N2  = p (number  of  pattern  points) 

AA  = a 


STEP  = step  size  e 

CUE  = (Z^  + Z,  1"^' 
A L 

EF  = F(^  ) 

P 
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ERR  = e(equation  24) 

GOM  = If  (4)  ) I 
'op' 

jk(x  cos4<  +y  sin({)  ) 

PAT  = e P " P 

V = (excitation  voltage) 

XLOAD  = (XjM 

ZA  = port  impedance  matrix  Z^(equation  12). 

Minimum  allocations  are  given  by 
COMPLEX  CUL(N),  EF(N2),  PAT(N,N2)  ZA(N*N) 

DIMENSION  A(N*N),  D(N),  G0M(N2) , NS(N),  U (N*N) , XL (N) , XLOAD (N) 


Subroutine  ERROR  Is  called  by  subroutine  ROSIEM  to  determine  a 
new  value  for  e for  given  and  . ROSIEM  searches  In  N 

variables  (N-1  reactive  loads  (the  driven  element  is  not  reactlvely 
loaded)  and  one  excitation  voltage). 

The  statement  "IF  (Nl-1)  10,  10,  11"  sets  [D]  = [I]  and  determines 
f(X°)  if  N1  = 1. 

DO  loop  14  is  the  main  search  loop  for  the  N variables,  Nested  DO 

loop  20  performs  two  functions.  The  first  is  to  increment  X^  and  the 

second  Is  to  determine  if  the  magnitude  of  the  current  load  value  has 

exceeded  500  ohms.  Once  the  load  value  has  reached  500  ohms.  It  has 

been  found  that  further  increases  in  this  load  value  toward  infinity 

will  not  noticeably  change  the  value  of  |F(({))|.  Therefore,  DO  loop 

20  will  determine  when  a load  magnitude  has  exceeded  500  ohms  and  sets 

the  load  magnitude  equal  to  A1  ohms  and  a sign  corresponding  to  the 

.2 

negative  of  the  former  reactance.  For  instance,  if  X^  = -550  ohms, 

2 

X becomes  + A1  ohms.  In  addition,  DO  loop  20  resets  the  value  of  E 
''  , , . ,N1 

(current  magnitude  of  e)  to  that  of  STEP  and  resets  N1  to  1 since 

will  be  incorrect  for  the  new  load  value.  Therefore  everytlme  a load 


value  goes  through  Infinity,  [D 


= [D^  = [I]. 
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Nested  statement  "GO  TO  (21 ,22,23,24) ,L"  corresponds  to  the  success 
and  failure  Increment  logic  for  each  search  cycle.  The  following  is  a 
summary  of  the  logic; 

a)  Statement  21  is  the  initial  screening  logic.  If  f(X^)  < f (X*^) 

/ . . N new  old  , „„ 

(the  first  step  is  a success),  e = ae  and  statement  23 

takes  over  for  subsequent  steps.  If  f(X^)  > f(X^)(the  first 

step  is  a failure),  and  statement  22  takes  over 

for  the  next  step, 

b)  Statement  22  is  used  if  the  first  step  was  a failure.  If 

2 1 

f (X  ) < f (X  ) (the  first  step  was  a failure  but  the  second  step 

was  a success) , = ae°^"^  and  statement  23  takes  over  for 

2 1 

subsequent  steps.  If  f (X  ) > f (X  ) (first  two  steps  were 
failures)  , = -ge°^^  and  statement  24  takes  over  for  sub- 

sequent steps. 

c)  Statement  23  is  used  if  the  first  step  was  successful  or  if 

the  first  step  was  a failure  and  the  second  step  was  success- 
ful. If  f(X^)  < f(X^“^),  = ae°^^.  If  f(X^)  > f(X^"^), 

nGw  old 

e = e /a  and  the  search  is  ended  for  this  variable  for  the 
present  search  cycle.  In  other  words,  statement  23  continues 
to  take  steps  as  long  as  they  are  successful  until  the  first 
failure  is  recorded. 

d)  Statement  24  is  used  if  the  first  two  steps  are  failures. 

If  f (xS  f (X*^) , the  search  is  ended  for  this  variable. 

If  f(X^  > f(X^)  (continuing  failures),  = -6e”^^  (a 

maximum  of  20  (-B)  steps  can  be  taken  before  the  search  ends 
in  this  variable). 

The  variable  NS(M)  records  the  number  of  steps  taken  for  variable 
number  M.  If  NS(M)  is  positive,  NS(M)  will  indicate  the  number  of  steps 
taken  in  the  successful  direction.  If  NS(M)  is  negative,  NS(M)  will  in- 
dicate the  number  of  steps  taken  in  the  failure  direction.  If  NS(M)  equals 
-20,  statement  31  will  discontinue  the  search  in  this  variable  until  a new 

direction  vector  d is  determined. 

1 

The  variable  D(M)  is  the  final  step  magnitude  for  the  current  search 


cycle. 


i 
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DO  loops  40  and  42  are  the  main  loops  for  finding  by  the 

Gram-Schmldt  orthogonallzatlon  procedure. 


3) 


Listing  of  Subroutine  ROSIEM. 
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SUBPrUTINf  ROSI  EM(N,N1  ,N2  ,AA,  STtPtCUL  , [ = ,GCM,PAT  ,V  ,XL  ■'Ar  , ZA 

1) 

C3MPLEX  CUL (7 ) ,EF(36) ,PAT(7,36), ZA(A9) 

D I ME  NS  I ON  A (A9)  ,0( 7) ,GOM(3b)  , NS ( 7 ) , U ( 49 ) , XL ( 7 ) ,XLCAD( 7) 

N3=N-1 

A1=30C. 

IF(Ni-l)  10,10,11 

10  J1=0 

DP  12  J=1 ,N 
J2=J1*J 
DP  13  1=1, N 
J1=J1+1 
13  U( J1  )=P. 

12  U(J2)=1. 

CALL  FRRPP  ( N,N2,FF  ,F  RP. , GpM.p/iT  , V.XLOAD,  ZA  ) 

11  DO  14  H=1 ,N 
NS(M) =0 

L=  1 

E = STf  P 
LRR1=ERP 
15  J1=M 
K = 2 

DP  20  1=1, N3 
XL(K)=XLOAP  (K  )+r,.u(  J1  ) 

IF(APS(XL{K)).GT,501.)  N1=0 
IF(AeS(XL(K)),GT.5ri.)  F=SIGN(STEP,XL(K)) 

IF (ABSCXL IK) ) .GT.501, ) XL  CAD ( K ) = -S I Gh ( A 1 , XL ( K ) ) 

] F<ABS(XL  (K  ) ) .GT.501.  ) XL  ( K ) =- S I GF  ( A 1 ,X  L ( K ) ) 

K=K  + 1 

20  Jl=Jl-*-N 
XL ( 1 )=0. 

V= V+F“U(J1 ) 

ERR2=ERR 

CALL  FRRUR(N,N2,EF,ERR,GnM,PA'",V,XL  ,ZA  ) 

DEPR=EPR-ERR2 

GO  TO  ( 21 ,22, 23,24) ,L 

21  L=3 

IF(DEPP)  25,25,26 
26  L = 2 

ERP=FPR2 
C=-STFP 
GO  TP  15 

22  L = 3 

IF(OEFP)  25,25,30 
3C  L = 4 

31  IF(NS(M)+20)  32,32,33 

33  E=-BB*F 

NS (M)=NS(M)-1 
GO  TP  15 
25  F=AA»E 

NS ( M) =NS (M) ♦! 

GO  TP  15 

23  IF(DFFR)  25,25,34 

34  E=F/AA 
FPP=FRP2 
GO  TO  32 

24  DFPR=FRR1-FRR 
IF(DFPP)  31,31,32 


32  J1=M 
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DT  35  1=1, N3 

XLrAD(K)  = XLn4D(K)  + 6*U(  Jl) 
K = K + 1 

35  J1=J1+N 

V = V + F=‘-U<  Jl) 

14  D(M)=F 

Jl  = N’*N-*-l 
DO  40  I=1,N 
J1=J1-1 

A( Jl) =D(N)»U( Jl) 

GO  TO  (4C I ,N 
J3  = N 

DO  41  M=2,N 
J2=J1 
J1=J1-1 
J3=J3-1 

41  A(  J1)=A(J2)+0(J3)''U(J1) 

40  CONTINUE 

DC  42  M=: ,N 
J1=M 

DO  43  I=1,N3 
XL(I  )=A(J1) 

43  J1=J1+N 

V = A( J]  ) 

GO  TO  (44),M 

MM=M-1 

DO  45  J=1,MM 

Z1=0. 

J1=J 
J2  = M 

DO  50  K=1,N 
Zl  = Zi-*-A(  J2)>^U(  Jl) 

J1=J1+N 

50  J2=J2+N 
J1  = J 
K = 2 

DO  51  1=1, N3 
XUK)=XL(K)-ZI*U(J1  ) 

K = K+1 

51  J1=J1+N 
V=V-Z1*U( Jl ) 

45  continue 

44  C=0. 

K = 2 

DO  52  1=1, N3 
C=C*XL(K)*XL( K ) 

52  K=K+1 
C=C-i-V>  V 
C=  l./SQFT (C ) 

Jl=r^ 

K = 2 

DO  53  1 =l ,N3 
U( J1)=C*XL(K) 

K=K+1 

53  J1=J1»N 
U( Jl  )=C^V 

42  CONTINuF 
P ETUPN 
FND 

Li 
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